Time-energy and time-entropy uncertainty relations in dissipative quantum dynamics 
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We derive exact relations and general inequalities that extend the usual time-energy uncertainty 
relations from the domain of unitary Hamiltonian dynamics to that of dissipative dynamics as de- 
scribed by a broad class of linear and nonlinear evolution equations for the density operator. We 
restrict our attention to intrinsic characteristic times of dynamical variables associated with the 
linear functionals Tr(p_F) of the density operator, as well as with the nonlinear entropy functional 
— fc B Tr(plnp). For non-dissipative dynamics, by using the Schroedinger inequality instead of the 
Heisenberg-Robertson inequality, we obtain a general exact time-energy uncertainty relation which 
is sharper than the usual Mandelstam-Tamm-Messiah relation tfAh > h/2. For simultaneous uni- 
tary/dissipative dynamics, the usual time-energy uncertainty relation is replaced by a less restrictive 
relation that depends on the characteristic time of dissipation, r, and the uncertainty associated 
with the generalized nonequilibrium Massieu-function operator which defines the structure of the 
dissipative part of the assumed class of evolution equations. Within the steepest-entropy-ascent 
dissipative quantum dynamics of an isolated system introduced earlier by this author, we obtain 
the interesting time-energy and time-entropy uncertainty relation (2tfAjj//i) 2 + (r^As /fe B r) 2 > 1. 
We illustrate this result and various other inequalities by means of numerical simulations. 

PACS numbers: 03.65.Ta,11.10.Lm,42.50.Lc,05.45.-a 



I. INTRODUCTION 



The time-energy uncertainty relation has remained an 
open and at times controversial issue throughout the his- 
tory of quantum theory. Several reviews are available 
on the pioneering discussions and the more recent devel- 
opments [ij. We are motivated by the recent revival of 
foundational questions and the development of dynami- 
cal theories that seek to build into the laws of quantum 
mechanics those of equilibrium and non-equilibrium ther- 
modynamics. Such revival is currently paralleled by the 
steady advancement of experimental techniques dealing 
with single ion traps @, qubits 0, neutron interferom- 
etry and a growing number of other developments 
all pointing at microscopic few-particle setups that nev- 
ertheless exhibit non-unitary dissipative dynamical be- 
haviour and call for the investigation of the role of gen- 
eral thermodynamic principles at the microscopic quan- 
tum level. Here we address the question of how the usual 
time-energy uncertainty relation, as interpreted accord- 
ing to the Mandelstam-Tamm-Messiah intrinsic-time ap- 
proach [j| based on unitary Hamiltonian dynamics, is 
modified by the presence of dissipation in the dynamical 
law or model. 



Let 7i (dimTL < oo) be the Hilbert space and H the 
Hamiltonian operator associated with a system in stan- 
dard Quantum Mechanics. We assume that the quantum 
states are one-to-one with the linear hermitian operators 
p on H. with Tr(p) = 1 and p > p 2 , and we assume a 



dynamical equation of the form 



dp 
dt 



P E{p)+E\p)p, 



(1) 



where E(p) is an operator-valued function of p that we 
may call the "evolution" operator and in general is non- 
hermitian. Without loss of generality, we write E = E+ + 
iE- where E + = (E + E*)/2 and E- = (E - E*)/2i are 
hermitian operators that, for convenience, we rename as 
AM(p)/2k B T(p) and H/h, respectively, so that Eq. (JTJ 
takes the form 



dp 
dt 



1 



2fc B r(p) 



{AM(p),p} , 



(2) 
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where [ • , • ] and { • , ■ } are the usual commutator and 
anticommutator, H (assumed independent of p, but not 
necessarily independent of time t) is identified with the 
Hamiltonian operator, h the reduced Planck constant, 
k B the Boltzmann constant; moreover, AAf(p) is a her- 
mitian operator-valued, possibly nonlinear function of p 
which, together with the positive definite, possibly non- 
linear functional r{p) of p, describes the dissipative dy- 
namics of the system, and is such that Tr[pAAf (p)] = 
as required to preserve p unit trace. 

The reason for considering a dynamical law of the 
form © is that the explicit expression of AM(p) that 
generates steepest-entropy-ascent (maximal entropy gen- 
eration) conservative dynamics of an isolated system 
compatible with all thermodynamics requirements [6j is 
known p| . We use it in Sections IVII and IVIII as an il- 
lustration of our general results, within a model for irre- 
versible relaxation in a Boltzmann gas of iV-level atoms. 
The extension to structured composite systems, which 
is nontrivial in view of the assumed nonlinearity of the 
dynamical law with respect to p, will be discussed else- 
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where. Additional discussion about the form of equation 
(J2J) is given in the Appendix. 

The steepest-entropy-ascent form of the operator func- 
tion AM(p) is discussed in Section IVll after introducing 
the necessary notation, in terms of the operator 



M(p) = S(p) 



9{p) 6{p) 



(3) 



where S(p) is the entropy operator [the precise defini- 
tion is given in Section ITTT1 Eq. (|33|l . and the function- 
al 9(p) and fi(p) are defined in Section IVll Eqs. H(i3|) - 
(|66|l ]. We call operator M(p) the generalized nonequilib- 
rium M as sieu- function operator, because at thermody- 
namic equilibrium its mean value belongs to the family of 
entropic characteristic functions introduced by Massieu 



i.e.. 



(M) c = (S) c 



(H)c , MN) C 



T 



T 



(4) 



where (S) c , (H) e , (N) c , T and fi are the equilibrium 
entropy, energy, amounts of constituents, temperature 
and chemical potentials, respectively. 

The specific physical interpretations of the uncertainty 
relations that we derive from dynamical law (J2J depend 
on the theoretical or modeling context in which it is as- 
sumed. In this article we limit the discussion to gener- 
alities in Sections [HI II I II IIVI Ivl and to illustrative con- 
siderations and numerical results valid within the sim- 
plest framework of stecpcst-cntropy-ascent conservative 
dynamics in Sections fVll and IVTT1 



II. GENERAL UNCERTAINTY RELATIONS 

We consider the space £(Tt) of linear operators on 7i 
equipped with the real scalar product 



(F\G) = Tr(FtG + G*F)/2 = (G\F) , 



(5) 



and the real antisymmetric bilinear form 

(F\G) = iTr(F^G-G^F)/2 = -(G\F) = (F\iG) , (6) 

so that for any (time-independent) hermitian F in C(Tl) 
the corresponding mean-value state functional can be 
written as (F) = Tr(pF) = Tt^Fy/p) = (y/p\y/pF), 
and can therefore be viewed as a functional of y/p, the 
square-root density operator, obtained from the spectral 
expansion of p by substituting its eigenvalues with their 
positive square roots. When p evolves according to Eq. 
JU, the rate of change of Tr(pF) can be written as 

dTr( P F)/di = Tr(F dp/dt) = 2 (y/pF\ y/pE(p)) . (7) 

In particular, for the evolution equation to be well 
defined, the functional Tr(p/) where I is the identity on 
7i must remain equal to unity at all times; therefore, 



dTV(p/)/dt = 2 (y/pl\ ^pE(p)) 
implies the condition 



0, in view of Eq. © 



Ol VpAM(p)) = 



(8) 



For F and G hermitian in C(H), we introduce the fol- 
lowing shorthand notation 



AF 

0~FG 

A F 

VFG 



F - Tr(pF)I , 

(AFAG) = (VpAF|V^AG) 



1 



Tr(p{AF, AG}) = a GF 



yfoFF = \J (AFAF) , 
({F,G]/2 i ) = (^pAF\^pAG) 
1 



^Tr(p[F,G\) 



Vfg 



-TjGF 



(9) 

(10) 
(11) 

(12) 



For example, we may write the rate of change of the mean 
value of a time-independent observable F as 



dTr(pF) ([F,H}/2i) (AFAM) 



dt 



h/2 



Vfh 
h/2 



Q~FAI 



(13) 

from which we see that not all operators F that com- 
mute with H correspond to constants of the motion, but 
only those for which (AFAM) = 0, i.e., such that yfpAF 
is orthogonal to both iyfpAH and y/pAM, in the sense 
of scalar product JSJ. For an isolated system, conserva- 
tion of the mean energy functional Ti(pH) requires an 
operator function AM(p) that maintains yfpAM always 
orthogonal to y/pAH, so that (AH AM) = for every p. 

From Schwarz inequality, we readily verify [!j the fol- 
lowing generalized Schrodinger uncertainty relation |lCj 



(AFAF) (AG AG) > (AFAG) 2 + ([F, G]/2i) 2 , (14) 



usually written in the form yf offo~gg — o~fg — \vfg\- Re- 
lation HI 4(1 obviously entails the less precise and less sym- 
metric Heisenberg-Robertson uncertainty relation 



(AFAF) (AG AG) > ([F,G]/2i 



(15) 



usually written in the form ApAc > |f?M?|- 

For further compactness, we introduce the notation 



tfg = °~fg / \/~off~&gg , 

CFG = VFG I \Jo-FFO-GG , 



(16) 



where clearly, rpG represents the cosine of the angle 
between the 'vectors' yfpAF and yfpAG in C(H), and 



' FG 



< 1. Inequality 1(14(1 may thus be rewritten as 



r FG + C FG — 1 



and clearly implies 



c fg — 



i + (4g/4g) 



< i - 4g < i 



(17) 



(18) 
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Next, for any hermitian F we define the characteristic 
time of change of the corresponding property defined by 
the mean value of the linear functional (F) = Tr (pF) as 
follows 



r F (p) = A F /\d(F)/dt\ 



(19) 



As is well known [1113; tf represents the time required for 
the statistical distribution of measurements of observable 
F to be appreciably modified, i.e., for the mean value (F) 
to change by an amount equal to the width Aj? of the 
distribution. 

Now, defining the nonnegative, dimensionless func- 
tional 



hA M 2k B rA H 



we rewrite ((13(1 in the form 

d(F)/dt = 2A F A H {c FH + a T r FM )/h 



(20) 



(21) 



and, substituting into (|19|l . we obtain the general exact 
uncertainty relation 



h/2 

TpA 



= \cfh + a T r FM \ 



H 



(22) 



For non-dissipative dynamics, a T - 
time-energy uncertainty relations 



0, Eq. Jill) yields the 



fe 2 /4 



l FH 



< 



1 



1 + (r FH /4 H ) 



< 1 



1 FH 



< 1 



(23) 



which entail but are more precise than the usual 
time-energy uncertainty relation, in the same sense as 
Schrodinger's relation (|14(l entails but is more precise 
than Heisenberg's relation 115|) . According to (|19fl . the 
last inequality in 1)23(1 implies that property (F) cannot 
change at rates faster than 2A F Ah /h. 

For dissipative dynamics let us first consider an observ- 
able A that commutes with H, so that ([A, H]/2i) = 
while (AAAH) ^ 0; in other words, an observable con- 
served by the Hamiltonian term in the dynamical law J5J , 
but not conserved by the dissipative term. Then Eq. 122(1 
yields the equivalent time-energy uncertainty relations 



h/2 

t a A 



a T \tam\ < a T 



H 



k B T(p) 

T A A 



\tam\ < 1 



M 



(24) 



(25) 



We note that while < 1, the value of a T depends 
on how r(p) is defined and, a priori, could well be larger 
than unity, in which case there could be some observables 
A for which taA# < h/2. If instead we impose that the 
functional r(p) be defined in such a way that a T < 1, i.e., 



than we obtain that even in dissipative dynamics the 
usual time-energy uncertainty relations are never vio- 
lated by observables A commuting with H . 

However, in general, if the dynamics is dissipative (r ^ 
oo ) there are density operators for which \c F n+a T tfm\ > 
1 so that tfAh takes a value less than h/2 and thus the 
usual time-energy uncertainty relation is violated. The 
sharpest general time-energy uncertainty relation always 
satisfied when both Hamiltonian and dissipative dynam- 
ics are active is (proof in Section llVjl 



h 2 /4 2 

< 1 + a + 2a T CMH 



T p&HH 

which may also take the equivalent form 

t f (thh t f omm t f A m A h cmh 
ft 2 /4 + klrHp) + k B T(p)h/4: 



> 1 



(27) 



(28) 



The upper bound in the rate of change of property (F) 
becomes 



A, 



&MM 



H 2 /4 + klr(p) + k B T(p) h/A 



(29) 



As anticipated, because the dissipative term in Eq. @ 
implies an additional dynamical mechanism, this bound 
((29(1 . valid for the particular nonunitary dynamics we are 
considering, is higher than the standard bound valid in 
unitary hamiltonian dynamics, given by 2ApAu/h. For 
observables commuting with H, however, ((25(1 provides 
the sharper general bound Aj?Am/';b''', solely due to dis- 
sipative dynamics, which is lower than 1(29(1 . 

Because in general \cmh\ < 1, l|2*S|l obviously implies 
the less precise relation 



ft 2 /4 

T F (THH 



< (l + «r) 2 



(30) 



However, as for the dynamics we discuss in Section IVT1 if 
the Massieu operator function AM(p) is a linear combi- 
nation (with coefficients that may depend nonlinearly on 
p) of operators that commute with either p or H , then it 
is easy to show that cmh — 0. Therefore, in such impor- 
tant case, ((28(1 becomes 



h 2 /A 

TpCTHH 



< 1 



(31) 



clearly sharper than 1(30(1 . If in addition r(p) satisfies 
(gSJl, then lEU implies t f A h > h/2^/2. 



III. RATE OF ENTROPY CHANGE 
CHARACTERISTIC TIME 



t( P ) > hA M /2k B A H 



We now consider the entropy functional (S) = 
( 26 ) Tr(pS) = -k B Tr(p\np) = -k B (^p\^pln(^p) 2 ) and its 
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rate of change, which using Eqs. J5J and (JSJ may be writ- 
ten as 

dTr( P 5)/d< = 2{^S\^E{p)) = (ASAM)/k B r 

= A S A M r SM /k B T (32) 
where, for convenience, we define the entropy operator 



S = -k B P p>0 \np 



(33) 



where P p >o is the projection operator onto the range of p 
. Interestingly, the rate of entropy change, being pro- 
portional to the correlation coefficient between entropy 
measurements and M measurements, under the assump- 
tions made so far, may be positive or negative, depending 
on how M is defined, i.e., depending on the specifics of 
the physical model in which Eq. © is adopted. 

The characteristic time of change of the entropy func- 
tional, defined as 



r s {p) = A s /|d(5)/di| 



(34) 



gives rise to the following equivalent exact time-energy 
uncertainty relations 



h/2 

t s A 



a T Wsm\ < a 7 



H 



k B T(p) 

T S A 



= \tsm\ < 1 



A I 



(35) 



(36) 



where tsm is defined as in (|I6J) using operators AM(p) 
and AS = S — (S). The physical interpretation of (|36|) is 
that the entropy cannot change in time at a rate faster 
than AsAm /k-aT, as immediately obvious also from (|32H . 

We notice from 135fl that if the dissipation time func- 
tional r(p) satisfies condition (|26[) then a T < 1 and, there- 
fore, the entropy change characteristic time ts satisfies 
the usual uncertainty relation tsAh > h/2 and the rate 
of entropy change cannot exceed 2AsAn/h. 

We conclude this Section by noting that, in general, 
the equality in l|35[) may be used to rewrite Relation 127|) 
in the form 



Or , , . V 1 + °r + 2a T C M H , , q _, 

\rsM\Ts < tf — ^ t f > (3') 



1 + a T 



1 + a r 



where the last inequality follows from | cmh | < 1 • This re- 
lation shows, on one hand, that the entropy change char- 
acteristic time ts is not necessarily the shortest among 
the characteristic times Tp associated with observables 
of the type (F) = Ti-(pF) according to the Mandelstam- 
Tamm definition (|19l) . On the other hand, it also shows 
that the left-hand side defines a characteristic-time func- 
tional 



T~UD 



1 



(38) 



which constitutes a general lower bound for all t^'s, and 
may therefore be considered the shortest characteristic 
time of simultaneous unitary /dissipative dynamics as de- 
scribed by Equation J3J). This observation prompts the 
discussion in the next Section. 



IV. SHORTEST CHARACTERISTIC TIMES 
FOR PURELY-UNITARY AND 
PURELY-DISSIPATIVE DYNAMICS 

The Mandelstam-Tamm definition l|19() of characteris- 
tic times has been criticized for various reasons (see for 
example Refs. 0> HE US) mainly related to the fact 
that depending on which observable F is investigated, as 
seen by inspecting (|23|l . the bound Tp > h/2 Ah may be 
very poor whenever c FH is much smaller than 1. 

Therefore, different attempts have been made to de- 
fine characteristic times that (1) refer to the quantum 
system as a whole rather than to some particular observ- 
able, and (2) bound all the particular r^'s from below. 
Notable examples are the characteristic times tes and 
t~lk, respectively defined by Eberly and Singh \'2 \\ and 
Leubner and Kiener |25| as follows 

The following definitions are based on the observa- 
tion that Ap may be interpreted as the norm of ^fpAF 
(viewed as a vector in C{TL)) in the sense that it equals 
•v/ (y/pAF\y/pAF) , therefore, we may use it to define the 
(generally non hcrmitian) unit norm vector in C(H) 

F p = y/pAF/ A F . (39) 

As a result, Eq. (|I3(I may be rewritten in the form 



where for shorthand we define the operator 



(40) 



C = 



. Ah&p A M M P 
1 h/2 k B r 



= 2^-pE(p) , (41) 



directly related [see Eq. Q] with the evolution operator 
function E(p) defined in the Introduction, which deter- 
mines the rates of change of all linear functionals of the 
state operator p, i.e., all observables of the linear type 
Tr(pF), by its projection onto the respective directions 

F P - 

Each characteristic time Tp can now be written as 

T F = A F /\d(F)/dt\ = l/\(F p \C)\. (42) 

Because F p is unit norm, | (F p \C) \ is bounded by the value 
attained for an operator F p that has the same 'direction' 
in £(H) as operator C, i.e., for 



F p - ±C/^{C\C) , 



(43) 



in which case \{F P \C)\ = y/(jO\Cj = y/Tr{CW). Thus 
we conclude that, for any, F, 



i/V(c\cj< 



T F 



(44) 



and, therefore, we introduce the shortest characteristic 
time for the combined unitary-dissipative dynamics de- 
scribed by Eq. @, 



tud 



(45) 
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which bounds from below all rp's. From (|41|l and 
103), and the identities (iH p \iH p ) = (M p \M p ) = 1 and 
(i#p|Afp) = (M p \iH p ) = cmh we obtain 



< 



1 

'UD 



= (C|C) = 



"MM 



(1 



ft 2 /4 fcir 2 ^) 
- a?. + 2a T c M H) , 



fc B T(p) ft/4 
(46) 



which proves Relations l(2"7|l and 

For nondissipative (purely Hamiltonian, unitary) dy- 
namics the same reasoning (or substitution of t = oo, 
a T = in the above relations) leads to the definition of 
the shortest characteristic time 

t v = h/2A H , (47) 

with which the usual time-energy relation reduces to 

t f > t v . (48) 

Its physical meaning is that when the energy dispersion 
(or uncertainty or spread) Ah is small, tjj is large and Tp 
must be larger for all observables F, therefore, the mean 
values of all properties change slowly |27|. i.e., the state 
p has a long lifetime. Conversely, states with a small 
energy spread cannot change rapidly with time. States 
that change rapidly due to unitary dynamics, necessarily 
have a large energy spread. 

Another interesting extreme case obtains from Eq. (J2J 
when AM(p) is such that the condition [p, H] = implies 
[AM (p) , H] = for any p, as for the steepest-entropy- 
ascent dynamics discussed in Sections lVll and lVlTl In this 
case, it is easy to see that if the state operator p com- 
mutes with H at one instant of time then it commutes 
with H at all times and, therefore, the entire time evo- 
lution is purely dissipative. Then, the reasoning above 
leads to the definition of the shortest characteristic time 



td = k B r/A 



M 



(49) 



It s noteworthy that Tp can be viewed as the character- 
istic time associated not with the (generally nonlinear) 
Massieu functional (M) = Tr(pM(p)) but with the lin- 
ear functional (A) = Tr(pA) corresponding to the time- 
independent operator A which at time t happens to co- 
incide with M{p{t)). 

For purely dissipative dynamics, the bound Tp > 
td = k s r/ Am implies that when Am/^bT, i-e., the ra- 
tio between the spread in our generalized nonequilibrium 
Massieu function and the dissipation time functional, is 
small, then Tp> is large and Tp must be larger for all ob- 
servables F, therefore, the state p has a long lifetime. 
Conversely, if some observable changes rapidly, Tp is 
small and since Tp, must be smaller, we conclude that 
the spread Am [more precisely, the ratio AM{p)/k B T(p)] 
must be large. 



In terms of t\j and Tp, we can rewrite l|20|) . (|36|l and 
gSJl as 



1 

TS 



tu/t d 



FsmI < J_ 

TD 



TD 



(50) 



(51) 







TFM 




V TJJ 


TD , 


< 


1 


1 




T UD 




< 


(" + 




\TU 


TD J 



2 cmh 
tutd 



(52) 



Eq. (|51|l implies that the entropy cannot change rapidly 
with time if the ratio AM(p)/k B T(p) is not large. The 
first equality in l|52|l follows from (F p \iH p ) = cph and 
(Fp\M p ) = rpM, which also imply that Eq. (|40|l may 
take the form 



d(F) 



TFM 
TD 



dt \ Tp 

and operator C defined in (|41|l takes also the forms 



(53) 



C 



.H, 

TU 



Ml 
TD 



,_.^AH + ^AM^ 



A h tu A m td 



and its norm is \/l/r^ + l/r|, + 2c M h/tutd- 

Similarly, the rate of entropy change (|32|l takes the 
form 



d(S) 
dt 



— (S P \M p ) 
td 



Ag r SM 

TD 



(55) 



which, because |tsm| < 1, implies the bounds [equivalent 
to dnnj and 1(51)1] . 

-^<^<^. (56) 

TD dt TD 



V. OCCUPATION PROBABILITIES 

An important class of observables for a quantum sys- 
tem are those associated with the projection operators. 
For example, for pure states evolving unitarily [27j, the 
mean value (P) = TV(p(i)P) where P = \<f>o){<f>o\ = p{0) 
represents the survival probability of the initial state, and 
is related to several notions of lifetimes [27j • 

We do not restrict our attention to pure states, and 
we discuss first results that hold for any projector P 
associated with a yes/no type of measurement. Let 
P = P' = P 2 be an orthogonal projector onto the g- 
dimensional subspace PTC oiTL. Clearly, g = Tr(P), the 
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variance (APAP) = p{\-p) where p = (P) = Tr(pP) 
denotes the mean value and represents the probability 
in state p of obtaining a 'yes' result upon measuring the 
associated observable, and the characteristic time of the 
rate of change of this occupation probability is defined 
according to fT§|) by 



1 

TP 



= 2 



= 2 



\dp/dt\ 

arcsin( v ^) 



arccos(Vp) 



< 



(57) 



where the inequality follows from l|46|) . Therefore, 

Id , ^ 1 

-5— < ^7 arccos(V^) < — — , 
It ud <tt It jjd 

or, over any finite time interval of any time history p(t), 



(58) 



arccos 



(VpW)) 



arccos 



< 



2T UD (t') 



(59) 

This result generalizes the results on lifetimes obtained 
in pi} where the focus is restricted to full quantum decay 
\p(oo) w 0] of an initially fully populated state [p(0) ~ 1] 
and tjj (here tud) is assumed constant during the time 
interval. It is also directly related to some of the results 
in where a number of additional inequalities and 
bounds on lifetimes are obtained for unitary dynamics, 
and may be straightforwardly generalized to the class of 
simultaneous unitary/dissipative dynamics described by 
our Eq. ©. 

Because p (1 — p) attains its maximum value when p = 
1/2, we also have the inequality 



dp 

at 



< 



1 



2tjjd 



(60) 



which, analogously to what noted in |26j, implies that 
no full decay nor full population can occur within a time 
2tud> so that this time may be interpreted as a limit to 
the degree of instability of a quantum state. 

Next, we focus on the projectors onto the eigenspaces 
of the Hamiltonian operator H, assumed time- 
independent. Let us write its spectral expansion as 
CnPe^ where e n is the n-th eigenvalue and P^ 
the projector onto the corresponding eigenspace. Clearly, 
PPw = enPe^, P^Psm = <WA, , 9n = Tr(Pj is the de- 
generacy of eigenvalue e„, p n = (P ?l ) = Tr^P^J the 
occupation probability of energy level e n , (AP^AP^) = 
Pn ($nm — Pm) the covariance of pairs of occupations, and 
(AP^ AP^) = p n (1 — p n ) the variance or fluctuation of 
the n-th occupation. Because [P^,if] = 0, cp^ = 0, by 
(1551) we have 



dp n 
dt 



= A, 



TD 



and the corresponding characteristic time is 



1 



TD 



TD 



(61) 



(62) 



Energy level occupation probabilities p n are used in 
Section IV11I for numerical illustration/validation of in- 
equalities (ff)2*|) within the steepest entropy ascent dynam- 
ical model outlined in the next Section. 



VI. EXAMPLE. 
STEEPEST-ENTROPY-ASCENT, 
CONSERVATIVE DISSIPATIVE DYNAMICS 

So far we have not assumed an explicit form of AM (p) 
except for condition JSJl that maintains p unit trace. In 
this section, we illustrate the above results by further as- 
suming steepest-entropy-ascent, conservative dissipative 
dynamics as obtained by assuming for our generalized 
nonequilibrium Massieu operator function the expression 



AM(p) = AS - AH'(p)/6(p) , (63) 

where S is the entropy operator defined in (J2HJ), 

AH'(p) = AH- n(p) • AN , (64) 

H is the Hamiltonian, N = {Nx, . . . , N r } a (possibly 
empty) set of operators commuting with H, that we call 
non-Hamiltonian generators of the motion (for example, 
the number-of-particles operators or a subset of them, or 
the momentum component operators for a free particle), 
such that operators ^fpAH and y/pAN are linearly in- 
dependent, 0(p) and /x(p) = {pi(p), . . . ,p r (p)} a set of 
real functionals defined for each p by the solution of the 
following system of linear equations 

r 

(AS AH) + (AN.AH) p t = (AH AH) , (65) 
i=i 

r 

(ASANj) 6 + (ANiANj) = (AHANj) ,(66) 

i=l 

which warrant the conditions that (AH AM) — and 
(ANjAM) = 0, and hence that the mean values Tr(pH) 
and Tr(pN) are maintained time invariant by the dissi- 
pative term of the equation of the motion. 

Operators y/pAH' and yfpAM are always orthogonal, 
in the sense that (AMAH 1 ) = for every p. It follows 
that, in general, (AS AH') = (AH'AH')/9, 



(AS AM) = (AM AM) = (AS AS) - 



(AH' AH') 
2 (P) 



>0 , 
(67) 

and hence the rate of entropy generation ll-'il'll is always 
strictly positive except for (AM AM) = (which occurs 
iff yfpAM = 0), i.e., for v /p^'AS , n d = (y r p~2AH - fi nd ■ 
^/p n dAN)/# n d, for some real scalars 8 n d and /^ n( j, that 
is, for density operators (that we call non-dissipative) of 
the form 



Pexp[-(g-Mnd-N)/Mnd]£ 
Pnd TrPcxp[-(p-- Atnd -N)/fc B ' 1 



(68) 
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where B is any projection operator on 7i (B 2 
The nonlinear functional 



B) 



8(p) = 



&H'H' 



&sh< A s r S H< 



(69) 



may be interpreted in this framework as a natural gener- 
alization to noncquilibrium of the temperature, at least 
insofar as for t — ► +00, while the state operator p(t) 
approaches a non-dissipative operator of form Q68[l. 9(t) 
approaches smoothly the corresponding thermodynamic 
equilibrium (or partial equilibrium) temperature # nt j- 

Because here H always commutes with M , cmh = 
and (M\iH) = 0, which means that ^fpAM(p) is always 
orthogonal to i^/pAH. This reflects the fact that on the 
entropy surface the direction of steepest entropy ascent 
is orthogonal to the (constant entropy) orbits that char- 
acterize purely Hamiltonian (unitary) motion (in which 
the entropy is maintained constant by keeping invariant 
each eigenvalue of p). 

Inequality i)67|) . which follows from r| A/ < 1, implies 
that (Jmm < oss and < tsm = Am / 'As < 1 or, equiva- 
lently, 



t k = k B T/A s < t d 



(70) 



where for convenience we define the characteristic time 
tk, which is simply related to the entropy spread, but 
cannot be attained by any rate of change, being shorter 
than T£>. In addition, we have the identities 



r SM 



VMM 



K 



D 



TS 



1 



1 



'SH' i 



and, from rg H , < 1, the bounds 



1*1 > 



A 



H' 



< I < As 



e - a 



H' 



(71) 



(72) 



where the equality \9\ — Ah> /As holds when and only 
when the state is non-dissipative [Eq. jSBjl]. Additional 
bounds on our generalized nonequilibrium temperature 9 
obtain by combining l|71|) with the inequality 4rg M (l — 
t Sm) — 1 (which clearly holds because r| M < 1), to obtain 
^ r SM r SH< — 1 anc ^' therefore, 



2A M A H , 



\6\ass 



< 1 or 



ess 



1 

< - < 



oss 



2A M A W ~ ~ 2A M A H , 

(73) 

At equilibrium, Am = and l|73[l implies no actual 
bound on 9, but in nonequilibrium states bounds l|73|) 
may be tighter than (|72|l , as illustrated by the numerical 
example in Section IVlII 

Notice that whereas in steepest entropy ascent dynam- 
ics tk'is always shorter than td and obeys the identity 



tstk 



ID ) 



(74) 



in general it is not necessarily shorter than td and obeys 
the identity 



Am t d 

As TSTK 



\rsM\ 



(75) 



In summary, we conclude that within steepest-entropy- 
ascent, conservative dissipative quantum dynamics, the 
general uncertainty relations l|28|l , i|35(l and H36JI that con- 
stitute the main results of this paper, yield the time- 
energy /time-Massieu uncertainty relation 



tfA h 
H/2 



M 



TpA 

k B T(p) 



> 1 or 4 + -4 > 1 , (76) 



1 D 



which implies the interesting time-energy/time-entropy 
uncertainty relation 



T F A H \ 2 ( TpAg N " 



h/2 



k B T(p) 



> 1 or + ^ 1 - ( 77 ) 



•U K 



and the time-entropy uncertainty relation 

t k k B T{p) 2 

— = ~ a = r SM < r SM < 1 , 

ts TsAs 



(78) 



which implies that the rate of entropy generation never 
exceeds crss/k B T, i.e., 

d(S) _ _ k Axr(plnp) - ° MM < ^ sAm < ass 
dt dt k B r ~ k B T ~ k B T 

(79) 

If in addition the dynamics is purely dissipative, such as 
along a trajectory p(t) that commutes with H for every 
t, then G7} may be replaced by the time-entropy uncer- 
tainty relation 



tk 
tf 



k B T(p) 
t f As 



< 1 . 



(80) 



As shown in Refs. [3, the dissipative dynamics gener- 
ated by Eq. J2J) with AM (p) as just defined and a time- 
independent Hamiltonian H: (i) maintains p(t) > p 2 (t) 
at all times, both forward and backward in time for any 
initial density operator p(0) (see also 0]); ( n ) maintains 
the cardinality of p(t) invariant; (iii) entails that the en- 
tropy functional is an ^-function in the sense defined in 
[13J and therefore that maximal entropy density opera- 
tors obtained from l|68[) with B — I are the only equilib- 
rium states of the dynamics that are stable with respect 
to perturbations that do not alter the mean values of the 
energy and the other time invariants (if any): this theo- 
rem of the dynamics coincides with a well-known general 
statement of the second law of thermodynamics : (iv) 
entails Onsager reciprocity in the sense defined in [15| ; (v) 
can be derived from a variational principle |12|. equiva- 
lent to our steepest entropy ascent geometrical construc- 
tion , by maximizing the entropy generation rate sub- 
ject to the Tr(p), Ti-(pH), and Tr(pN) conservation con- 
straints and the additional constraint ( K ^fpE\ y fpE) = c(p) 

We finally note that assuming in Eq. J5J , in addition to 
AM(p) given by l|63[l . the nonlinear relaxation time r(p) 
given by (|26|l with strict equality, we obtain the most dis- 
sipative (maximal entropy generation rate 17]) dynamics 
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in which the entropic characteristic time ts [Eq. (f 341) ] is 
always compatible with the time-energy uncertainty re- 
lation tsAjj > h/2 and the rate of entropy generation is 
always given by 2 Am Ah /h. 

The physical meaning of relations l|35|). (|5uTl . 

(|78|) . I|77ll are worth further investigations and experi- 
mental validation in specific contexts in which the dissi- 
pative behavior is correctly modeled by a dynamical law 
of form (0), possibly with AM(p) of form Ipfy. One such 
context may be the currently debated so-called "fluctu- 
ation theorems" 01 whereby fluctuations and, hence, 
uncertainties are measured on a microscopic system (op- 
tically trapped colloidal particle 0, electrical resistor 
|20|) driven at steady state (off thermodynamic equilib- 
rium) by means of a work interaction, while a heat in- 
teraction (with a bath) removes the entropy being gen- 
erated by irreversibility. Another such context may be 
that of pion-nucleus scattering, where available experi- 
mental data have recently allowed partial validation [2l| 
of "entropic" uncertainty relations 123. Yet another is 
within the model we propose in Ref. j2j] for the descrip- 
tion of the irreversible time evolution of a perturbed, iso- 
lated, physical system during relaxation toward thermo- 
dynamic equilibrium by spontaneous internal rearrange- 
ment of the occupation probabilities. We pursue this 
example in the next section. 



VII. NUMERICAL RESULTS FOR 
RELAXATION WITHIN A DILUTE 
BOLTZMANN GAS OF iV-LEVEL PARTICLES 

To illustrate the time dependence of the uncertainty 
relations derived in this paper, we consider an isolated, 
closed system composed of noninteracting identical par- 
ticles with single-particle eigenstates with energies for 
i = 1, 2,. . . , N, where N is assumed finite for simplicity 
and the e^'s are repeated in case of degeneracy, and we 
restrict our attention to the class of dilute-Boltzmann- 
gas states in which the particles are independently dis- 
tributed among the N (possibly degenerate) one-particle 
energy eigenstates. This model is introduced in Ref. 
|23|. where we assume an equation of form with 
AM(p) given by l|63l) with the further simplification that 
AH'(p) = AH so that our generalized noncquilibrium 
Massieu-function operator is 



M(p) = S- H/9(p) 



and, therefore, 



AM(p) = AS - AH/9{p) 



(81) 



(82) 



For simplicity and illustrative purposes, we focus on 
purely dissipative dynamics by considering a particular 
trajectory p(t) that commutes with H at all times t, as- 
suming that H is time independent and has a nondegen- 
erate spectrum. As a result, the energy-level occupation 
probabilities p n coincide with the eigenvalues of p, and 



the dynamical equation reduces to the simple form [2 

(S) e n - (H ' 





1 


dt 


T 


where 






(S) 




(H) 








Pn lnp„ 



Pn 



Pn 



-/c B ^Pnhip„ , 

n 

5 ' p« e « 1 
n 

J2p n el -(H) 2 , 

n 

cths = -k B ^p n e n lnp n - (H)(S) 



(83) 

(84) 

(85) 

(86) 
(87) 

(88) 



To obtain the plots in Figures 1-4, that illustrate the 
main inequalities derived in this paper for a sample tra- 
jectory, we consider an initial state with cardinality equal 
to 4, with nonzero occupation probabilities only for the 
four energy levels e\ — 0, e 2 = u/3, = u/3, and e 4 — u, 
and with mean energy (H) = 2u/5 (u is arbitrary, with 
units of energy). Moreover, as done in [2j|, we select an 
initial state p(0) at time t — such that the resulting tra- 
jectory p(t) passes in the neighborhood of the partially 
canonical nondissipative state p that has nonzero oc- 
cupation probabilities only for the three energy levels e\ , 
e 2 , and 64, and mean energy (H) = 2u/5 (p„ dl = 0.3725, 
p&a = 0.3412, <j 3 = 0, p£ d4 = 0.2863, 0* d = 3.796 u/k B 
). As shown in Figure 1, during the first part of the tra- 
jectory, this nondissipative state appears as an attractor, 
an approximate or 'false target' equilibrium state; when 
the trajectory gets close to this state, the evolution slows 
down, the entropy generation drops almost to zero and 
the value of gets very close (3.767 u/k B ) to that of 0^ d ; 
however eventually the small, but nonzero initial occupa- 
tion of level e% builds up and a new rapid rearrangement 
of the occupation probabilities takes place, and finally 
drives the system toward the maximal entropy state /?^ d 
with energy (H) = 2u/h and all four active levels occu- 
pied, with partially canonical distribution p^ dl = 0.3474, 
pl c d2 = 0.2722, = 0.2133, ? Z 4 = 0.1671, and char- 
acterized according to (|68|) by the (partial equilibrium) 
temperature 6^ d = 1.366 u/k B . 

The trajectory is computed by integrating Eq. I|83|) nu- 
merically, both forward and backward in time, starting 
from the chosen initial state p(0), and assuming for Fig- 
ures 1 and 2 that the dissipation time r is a constant, 
and for Figures 3 and 4 that it is given by l|26|) with 
strict equality (a T = 1, t d =tu), i.e, assuming 



h/2 A 



M 



<?SS 



k B A H 
1.2 \ " 



h/2 



Pn(lnp„ 



<?ss 

(sy 



1 



(89) 
(90) 



The system of ordinary differential equations (|83|l is 
highly nonlinear, especially when r is assumed accord- 
ing to (|89|) . nevertheless it is sufficiently well behaved to 
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allow simple integration by means of a standard Runge- 
Kutta numerical scheme. Of course, we check that at all 
times -co < t < oo each p n remains nonnegative, ^2 n p n 
remains equal to unity, ^2 n Pn^n remains constant at the 
value 2u/5 fixed by the selected initial state, and the rate 
of change of (S) is always nonnegative. 

In each Figure, the top subfigure shows for ease of com- 
parison the plots of the four nonzero occupation proba- 
bilities as functions of dimensionless time: t/r, in Figures 
1 and 2; ut/H, in Figures 3 and 4. The dots on the right 
represent the maximal entropy distribution, p n (+oo) = 
p^ d ; the dots at the left represent the lowest-entropy 
or 'primordial' distribution, p n (— oo) = p£ d , which for 
the particular trajectory selected here, corresponds to 
a nondissipative state p l * d that has only two occupied 
energy levels, e\ and e±, with probabilities p^ dl = 0.6 
and p^ di = 0.4 (and temperature p^ d = 2.466 u/k B ); in 
fact the system has no lower entropy states p that com- 
mute with H, have energy 2u/5, and have zero occupa- 
tion probabilities p n for every n > 4 |23j . The dots in 
the middle represent the nondissipative state p^ d which 
appears as the false target state during the first part of 
the trajectory, plotted at the instant in time when the 
entropy of the time- varying trajectory is equal to the en- 
tropy of this distribution. 

It is interesting to observe from Figures 1 (bottom sub- 
figure), and Figure 2 (second subfigure) that during the 
early part of the trajectory, td almost exactly coincides 
with Tp e2 while in the late part it almost exactly coin- 
cides with Tp , and the switch occurs when the trajec- 
tory slows down in the neighborhood of the 'false target' 
nondissipative state. 

In Figure 1, the second subfigure shows the time de- 
pendence of the dimensionless entropy (S)/k B ; the third 
subfigure shows its rate of change (proportional to omm) 
and compares it with ass and (Jhh/6 2 , to illustrate re- 
lation H67|l : the fourth shows the time dependence of 
our generalized 'nonequilibrium temperature' 9 (prop- 
erly nondimensionalized) and compares it with Ajy/Ag 
and 2AmAjj/oss to illustrate relations (f72*)) and l(73): 
the fifth subfigure shows the time dependence of 1/t d 
(which here is proportional to the square root of the rate 
of entropy generation, third subfigure) and compares it 
with 1 /tk and 1 /t$ to illustrate relations (|70|l and (15 1 1) ; 
the sixth subfigure shows 1/rp for each of the four oc- 
cupation probabilities and compares them with 1/td to 
illustrate relation l|62|) . which for this particular trajec- 
tory has the feature we just discussed. 

In Figure 2, the second subfigure illustrates again rela- 
tion (|62(1 for each of the four observables p n = (Pe^)', the 
third subfigure illustrates the time-entropy uncertainty 
relation l|8U|) for the same observables; the fourth illus- 
trates inequality (|60[) : the fifth illustrates relations i|51[) 
and (25 . 

Similar remarks hold for Figures 3 and 4, where it 
is noteworthy that most qualitative features remain the 
same, except for the almost singular behavior near canon- 
ically distributed nondissipative states, where Am ap- 



proaches zero and so does the dissipative time r assumed 
in this case according to (|%9*)l . The approach to (partial) 
equilibrium in this case is not exponential in time as for 
t = const. This puzzling behavior suggests that assump- 
tion l|89l) may hardly be physically sensible. However, as 
already noted after l|24|) , it represents an extreme behav- 
ior, i.e., the minimum dissipative time functional r(p) 
by which observables which commute with H, like the 
occupations P^ , never violate the usual time-energy un- 
certainty relations Tp Ah > h/2, even though their time 
dependence is not determined here by unitary dynamics 
but by purely dissipative dynamics. These usual time- 
energy uncertainty relations, Tp > tjj, are illustrated 
by the second subfigure of Fig. 4, because in this case 



VIII. CONCLUSIONS 

The Mandelstam-Tamm-Messiah time-energy uncer- 
tainty relation tfAh > h/2 provides a general lower 
bound to the characteristic times of change of all observ- 
ables of a quantum system that can be expressed as linear 
functionals of the density operator p. This has been used 
to obtain estimates of rates of change and lifetimes of 
unstable states, without solving explicitly the time de- 
pendent evolution equation of the system. It may also 
be used as a general consistency check in measurements 
of time dependent phenomena. In this respect, the exact 
relation and inequalities l|22|) [that we derive for standard 
unitary dynamics based on the generalized Schrodinger 
inequality l|14|) ] provide a more general and sharper chain 
of consistency checks than the usual time-energy uncer- 
tainty relation. 

The growing interest during the last several decades 
in quantum dynamical models of systems undergoing ir- 
reversible processes has been motivated by impressive 
technological advances in the manipulation of smaller 
and smaller systems, from the micrometer scale to the 
nanometer scale, and down to the single atom scale. The 
laws of thermodynamics, that fifty years ago were in- 
variably understood as pertaining only to macroscopic 
phenomena, have gradually earned more attention and 
a central rule in studies of mesoscopic phenomena first, 
and of microscopic phenomena more recently. In this pa- 
per we do not address the controversial issues currently 
under discussion about interpretational matters, nor do 
we attempt a reconstruction and review of the different 
views, detailed models and pioneering contributions that 
propelled this fascinating advance of thermodynamics to- 
wards the realm of few particle and single particle sys- 
tems. 

However, motivated by this context and background, 
we derive extensions of the usual time-energy uncertainty 
relations that extend their usefulness to studies of dissi- 
pative phenomena. We do this by focusing on a special 
but broad class of model evolution equations, that has 
been designed for the description of dissipative quantum 
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phenomena and for satisfying a set of compatibility con- 
ditions with general thermodynamic principles, and we 
derive in this framework various forms of time-energy and 
time-entropy uncertainty relations, and other interesting 
general inequalities, that should turn out to be useful as 
additional consistency check in measurements of time- 
dependent dissipative phenomena. Finally, we illustrate 
and discuss some of these relations with numerical re- 
sults obtained by integration of the nonlinear evolution 
equation introduced by this author for the description of 
steepest entropy ascent dynamics of an isolated system 
far from thermodynamic equilibrium. 



APPENDIX A: REASONS FOR NOT ASSUMING 
A KOSSAKOWSKI-LINDBLAD FORM OF THE 
EVOLUTION EQUATION 

With various motivations, fundamental or phenomcno- 
logical, a number of different generalizations of quan- 
tum dynamics have been proposed in which the evolu- 
tion equation for the density operator p does not con- 
serve the functional — Tr(plnp). In particular, in theo- 
ries of systems in contact with a heat bath, or subsystems 
of composite systems which as a whole evolve unitarily, 
a variety of successful model evolution equations for the 
reduced density operator have the Kossakowski-Lindblad 
[H form 

% = -l[H,p] + \ y E(2V}pV i -{V}V Jt p}) , (Al) 



where the Vj's are operators on Tt (each term within 
the summation, often written in the alternative form 
[Vj,pV^] + [VjP, Vj], is obviously traceless). Evolution 
equations of this form are linear in the density operator 
j^_and preserve its hermiticity, nonnegativity and trace 

For example, in a number of successful models of dissi- 
pative quantum dynamics of open subsystems, operators 
Vj arc in general interpreted as creation and annihila- 
tion, or transition operators. For example, by choosing 
Vj = c rs \r}(s\, where c rs are complex scalars and |s) 
eigenvectors of the Hamiltonian operator H , and defin- 



which, for the occupation probabilities p n = p nn , is the 
Pauli master equation 



ing the transition probabilities w rs 
(IA1I) becomes 



equation 



_(A2) 

or, equivalcntly, for the nm-th matrix element of p in the 
H representation, 

dPnm _ _i_ ( TP _ TT \ 

, . . Pnrn K-^n ) 

at n 

"F^ram ^ W nr p rr Pnm~Z ^ ' ( m rn "F w rm) i (A3) 



-TT = 2^ W nrPr - Pn ' 



(A4) 



In this paper, we consider a class of model evolution 
equations applicable not only to open systems but also 
to closed isolated systems, capable of describing, simul- 
taneously with the usual Hamiltonian unitary evolution, 
the natural tendency of any initial nonequilibrium state 
to relax towards canonical or partially-canonical thermo- 
dynamic equilibrium, i.e., capable of describing the irre- 
versible tendency to evolve towards the highest entropy 
state compatible with the instantaneous mean values of 
the energy, the other constants of the motion, and possi- 
bly other constraints. To avoid the severe restrictions im- 
posed by the linearity of the evolution equation, we open 
our attention to evolution equations nonlinear in the den- 
sity operator p. Therefore, it may at first appear natu- 
ral to maintain the Kossakowski-Lindblad form IjAlfl and 
simply assume that operators Vj are functions of p. This 
is true only in part for the evolution equation (J2J) that we 
assume. Indeed, our hermitian operator AM(p)/k B r(p) 
can always be written as — J2j Vj{p)Vj(p) and therefore 
our anticommutator term may be viewed as a generaliza- 
tion of the corresponding term in IjAlfl . 

However, in our equations and @ we suppress the 
term corresponding to ^2jVjpVj in <|A1|1 . The reason 
for this suppression is the following. Due to the terms 
Vj pVj, whenever the state operator p is singular, i.e., it 
has one or more zero eigenvalues, Eq. (|A1|) implies that 
these zero eigenvalues may change at a finite rate. This 
can be seen clearly from (|A4ll by which dp n /dt is finite 
whenever there is a nonzero transition probability w nr 
from some other populated level (p r ^ 0), regardless of 
whether p n is zero or not. When this occurs, for one 
instant in time the rate of entropy change is infinite, as 
seen clearly from the expression of the rate of entropy 
change implied by IjAlfl . 

^ = k^^Vjphip-VjpVjbxp) 

3 

= k B ^2(Vj)* nr (Vj) nr {p r - pn)\np r , (A5) 

jm 

where p r denotes the r-th eigenvalue of p and {Vj) nr the 
matrix elements of Vj in the p representation. 

We may argue that an infinite rate of entropy change 
can be tolerated, because it would last only for one in- 
stant in time. But the fact that zero eigenvalues of p 
in general could not survive, i.e., would not remain zero 
(or close to zero) for longer than one instant in time, is 
an unphysical feature, at least because it is in contrast 
with a wealth of successful models of physical systems 
in which great simplification is achieved by limiting our 
attention to a restricted subset of relevant eigenstates 
(forming a subspace of TL that we call the effective Hilbert 
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space of the system |6j). Such common practice mod- 
els yield extremely good results, that being reproducible, 
ought to be relatively robust with respect to including in 
the model other less relevant eigenstates. In fact, such 
added eigenstates, when initially unpopulated, are irrele- 
vant if they remain unpopulated (or very little populated) 
for long times, so that neglecting their existence intro- 
duces very little error. The terms Vj pVj, instead, would 
rapidly populate such irrelevant unpopulated eigenstates 
and void the validity of our so successful simple models, 
unless we deliberately overlook this instability problem 
by highly ad- hoc assumption, e.g., by forcing the V/s 
to be such that (Vj) nr — whenever either p n = or 
p r = 0, in which case, however, we can no longer claim 
true linearity with respect to p. 

To avoid the unphysical implications of this seldom 
recognized 0, 0] problem of linear evolution equations of 
form (|A1|) . we consider in this paper only equations of 
form (|2J|. We do not exclude that it may be interesting 
to investigate also the behavior of equations that include 
nonlinear terms of the form vj (p) p Vj (p) . However, at 
least when the system is strictly isolated, the operator- 
functions Vj(p) should be such that (Vj(p)) nr — when- 
ever either p n = or p r = 0. 

Another important general physical reason why we ex- 
clude terms that generate nonzero rates of change of zero 
eigenvalues of p, is that if such terms are construed so as 
to conserve positivity in forward time, in general they 
cannot maintain positivity in backward time. The view 
implicitly assumed when Eq. I|A1|I is adopted, is that 



the model is "mathematically irreversible" (a distinguish- 
ing feature if not a starting point of the theory of com- 
pletely positive linear dynamical semigroups on which it 
is based) , in the sense that neither uniqueness of solutions 
in forward time nor existence in backward time are re- 
quired (and granted). Such mathematical irreversibility 
of the initial value problem, is often accepted, presented 
and justified as a natural counterpart of physical irre- 
versibility. However, it is more related to the principle of 
causality than to physical irreversibility. The strongest 
form of the non-relativistic principle of causality, a key- 
stone of traditional physical thought, requires that future 
states of a system should unfold deterministically from 
initial states along smooth unique trajectories in state 
domain defined for all times (future as well as past). Ac- 
cepting mathematical irreversibility of the model dynam- 
ics, implies giving up such causality requirement. But it 
is not strictly necessary to describe physical irreversibil- 
ity, at least not if we are willing to give up linearity in- 
stead. The proof of this statement is our Eq. J5J) which, 
together with the additional assumptions made in Sec- 
tion ^T] to describe relaxation within an isolated system, 
is mathematically reversible, in the sense that it features 
existence and uniqueness of well-defined solutions both 
in forward and backward time, and yet it does describe 
physically irreversible time evolutions, in the sense that 
the physical property described by the entropy functional 
— fc B Tr(plnp) is a strictly increasing function of time for 
all states except the very restricted subset defined by Eq. 
(|68p. where it is time invariant. 
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FIG. 1: (Color online) See text; r assumed constant. 
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FIG. 3: (Color online) See text; r assumed according to 1891 . 



